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/ discuss the chemical and spectrophotometric evolution of galaxies over 
cosmological timescales and present a first attempt to treat both aspects 
in a chemically consistent way. In our evolutionary synthesis approach, 
we account for the increasing metallicity of successive generations of 
stars and use sets of stellar evolutionary tracks, stellar yields, spectra, 
etc. for various metallicities. This gives a more realistic description 
of nearby galaxies, which are observed to have broad stellar metalllicty 
y—i • distributions, as well as of young galaxies at high redshift. Selected re- 

sults are presented for the chemo- cosmological evolution of galaxies as 
compared to QSO absorption line observations and for the spectropho- 
tometric evolution of galaxies to very high redshifts. On cosmological 
timescales, intercations are important drivers of galaxy evolution. Ne- 
glecting dynamical aspects we study the effects of interaction-induced 
starbursts on the spectrophotometric and chemical evolution of galax- 
ies and briefly discuss the formation of star clusters and Tidal Dwarf 
Galaxies in this context. 

■ 

1 1 Introduction 

a: 

The evolution of galaxies on cosmological timescales has 3 aspects: the chem- 
ical, spectral, and dynamical evolution. In reality, all of them are intimately 
coupled as shown by the existence of several observational relations involving 
quantities from any two of these 3 domains. Examples are the Faber - Jack- 
son, Fundamental Plane, or luminosity - metallicity relations for ellipticals 
or the Tully - Fisher relation and trends of characteristic HII region abun- 
dances, average colours, emission line strengths, luminosities, and mass - to - 
light ratios with spiral type, i.e. with the transition from bulge-dominated to 
disk-dominated systems. 

While over short timescales or lookback times and, at least, for giant 
galaxies, the 3 aspects of galaxy evolution might be treated independently, any 
modelling of galaxy evolution over cosmological timescales definitely requires 
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a consistent treatment of all 3 aspects. Unfortunately, this is still too complex 
today, even for numerical simulations on powerful machines. 

Independently, all 3 aspects of galaxy evolution modelling have quite some 
tradition already. We are, however, still at the very beginning of a consistent 
understanding of galaxy formation and evolution over cosmological timescales 
including all 3 aspects. A consistent modelling of the chemical and dynamical 
evolution of galaxies is e.g. attempted by the group of G. Hensler in Kiel 
('Chemodynamics'). With the availability of stellar input physics for vari- 
ous metallicities, we developed in Gottingen a consistent description of the 
chemical and spectral evolution over cosmological timescales, i.e. coupled to a 
cosmological model. This chemically consistent evolutionary synthesis 
method and selected applications for the interpretation of galaxy observations 
over a wide redshift range will be presented here. 

While still being far from a consistent coupling of our chemically consis- 
tent chemical, spectrophotometric and cosmological model with a dynamical 
galaxy evolution model, we at least started to study the starbursts accom- 
panying strong interactions and mergers involving gas rich galaxies and their 
impact on the chemical and spectral evolution of galaxies. This brings us 
to some surprising features accompanying those strong bursts: the formation 
and evolution both of a new population of bright and massive young star clus- 
ters and of a new class of dwarf galaxies forming from recycled material along 
the tidal features of the merging giant galaxy pair (Tidal Dwarf Galaxies) . 

2 Chemically Consistent Evolutionary Synthe- 
sis 

2.1 Observational Evidence for Broad Metallicity Dis- 
tributions and Subsolar Metallicities 

Evidence for extended stellar metallicity distributions and, in particular, for 
the importance of subsolar abundance contributions has been accumulating 
over the last years. Stars in elliptical galaxies and bulges, e.g., show a metallic- 
ity distribution extending over more than a factor 10 from 0.06 < Z*/Z Q < 2.5 
with average metallicity slightly below solar (Mc William & Rich 1994), G 
dwarf stars in the solar neighbourhood cover a range from 0.15 < Z*/Z Q < 1.5 
with, again, the average metallicity being subsolar (Rocha-Pinto & Maciel 
1996). Characteristic HII region abundances (i.e. measured at 1 R e ) in spi- 
rals are observed to be in the range from ~ Z for early type spirals Sa/b 
through -C Z Q for late type spirals Sd (Oey et al. 1993, Zaritsky et al. 1994, 
Ferguson et al. 1998). Dwarf galaxies of all classes (dEs, dSphs, dls, BCDs) 
show - sometimes substantially - subsolar metallicities, both in their stellar 
populations and in their gas phase abundances (e.g. Richer & McCall 1995). 
So, while already in the local Universe, subsolar abundances are clearly pre- 
vailing on global galaxy scales - with the exception of some central regions 
of massive luminous galaxies - this is, of course, even more so the case in the 
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early Universe. 

From restframe UV stellar wind lines, redshifted into the optical for star 
forming galaxies at redshifts 3-4, stellar abundances around 1/10 Z© are 
derived (cf. Lowenthal et al. 1997, Trager et al. 1997). Damped Lya ab- 
sorbers - most probably the progenitors of present day galactic disks - fea- 
ture ISM abundances in the range 10~ 3 < Zism/Z© < 1 over a redshift range 
0.4 < z < 4.5 (e.g. Pettini et al. 1997, 1999, see also Sect. 5). 

These are only a few selected examples from a much longer and still in- 
creasing list of direct observational hints to the existence of significant metal- 
licity distributions in galaxies and to the importance of subsolar metallicities. 

2.2 Method and Input Physics 

From a very basic consideration it is immediately clear that any stellar system 
with star formation (SF) extending over more than the lifetime of the most 
massive stars (10 6 yr) is composite both in terms of age and metallicity. 

While star clusters are Simple Stellar Populations (SSPs) consisting of 
one stellar generation with one age and one metallicity, any galaxy is a 
Composite Stellar Population (CSP) with its stars spanning finite ranges 
both in age and metallicity. 

The Gottingen Evolutionary Synthesis Model starts from a gas cloud 
of given mass, assumes some Star Formation History (SFH) and Initial Mass 
Function (IMF) - the basic parameters of this kind of approach - and then, 
with some book-keeping algorithm, calculates the time evolution of the stellar 
population across the HRD from a set of stellar evolutionary tracks. With 
photometric calibrations for colours or absorption indices along all stellar evo- 
lutionary tracks the photometric evolution in terms of luminosities, colours, 
and absorption indices is obtained. Assigning stellar spectra to any point 
along the tracks and weighting them with the numbers of stars calculated at 
these points at any given time yields the time evolution of a synthetic galaxy 
spectrum. Integrating over the synthetic galaxy spectrum with response func- 
tion for any filter system (e.g. Johnson UBVRIJHK, HST F300W, F450W, 

F814W), luminosities in the respective passbands are obtained in the same 
way as for an observed galaxy spectrum. 

In this approach, evolutionary consistency is guaranteed and the time 
evolution can explicitly be studied. This opens the possibility to directly 
couple to a cosmological model and study the spectrophotometric evolution 
as a function of redshift over cosmological timescales (cf. F.-v.A. 1989). 

Evolutionary synthesis modelling per se accounts for the age distribution 
of the stars within a galaxy and its time evolution (for a review on synthesis 
methods see F.-v.A. 1994). To properly account for the metallicity distribu- 
tions in composite systems like galaxies, the abundance evolution in the gas 
has to be followed in order to describe successive generations of stars forming 
out of this gas with stellar evolutionary tracks, yields, and spectra appropriate 
for their respective initial metallicities, i.e. the gas metallicity at their birth. 
We call this approach chemically consistent (= cc). 
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A modified form of Tinsley's equations (Tinsley 1980) with stellar yields 
for SNII, SNI, PN, and stellar mass loss is solved to obtain the gas content G, 
the global metallicity Z of the ISM, as well as individual element abundances 
[Xj/H], and abundance ratios [X;/Xj]. We use sets of input physics for 5 
different stellar metallicities in the range 

-2.3 < [Fe/H] < +0.3 

with logarithmic element abundances with respect to solar defined by [X;/H] 
:= Zog(Xi/Xn) — log(Xf /X!^). These sets of input physics for the description 
of both the spectral and chemical evolution of model galaxies comprise stel- 
lar evolutionary tracks - covering all relevant evolutionary phases from ZAMS 
through the PN phase or SN explosion - lifetimes, and remnant masses from 
the Padova group (Bressan et al. 1993, Fagotto et al. 1994a, b, c). Evolu- 
tion of low mass stars is taken from Chabrier & Baraffe 1997. Stellar yields 
for a series of individual elements from 12 C through 56 Fe for massive stars 
(> 8 M Q ) are from Woosley & Weaver 1995 and yields for intermediate mass 
stars from van den Hoek & Groenewegen 1997. SNIa contributions to Fe, 
C, are included for the carbon deflagration white dwarf binary scenario as 
outlined by Matteucci & Greggio 1986. We caution that metallicity dependent 
stellar yields depend on -^r, explosion energies, remnant masses, etc. and the 
metallicity dependence of these factors is still poorly understood. SNIa yields 
are only available for Z Q (Nomoto et al. 1997, model W7). No important 
metallicity dependence is expected for SNIa yields except for a possible lower 
metallicity limit to the explosion (Kobayashi et al. 1998). Model atmosphere 
spectra, colour and absorption index calibrations are from Lejeune et al. 1997, 
1998, and Worthey et al. 1994, respectively. 

Although this cc approach goes significantly beyond what was possible 
before, it still is a simplification in the sense that all the input physics is only 
available for scaled solar abundance ratios. Abundance ratios in galaxies 
are determined by the metallicity dependent stellar yields, the SFH and IMF, 
and, in general, will be non-solar. I.e., stellar evolution and galaxy evolution 
themselves are implicitly coupled and the coupling depends on the SFH and 
IMF (cf. F.-v.A. 1998b). For a review on chemically consistent evolutionary 
synthesis see F.-v.A. 1999b. 

The spectrophotometric evolution of the stellar component of galaxies and 
the chemical evolution of ISM abundances both are studied not only as a 
function of time, but also - for any cosmological model as given by H , Q Q , A , 
and a redshift of galaxy formation zf - as a function of redshift. 

2.3 Model Parameters and Local Templates 

We use an IMF in the form given by Scalo 1986 with lower and upper 
mass limits of 0.08 and 85 M Q , respectively, and SFHs, \P(t), appropriate 
for the various spectral types of galaxies. For ellipticals we use a standard 
*(t) ~ c _t /**. For spiral types Sa ... Sc the Star Formation Rate (SFR) 
at any timestep is tied to the gas-to-total mass ratio, ^(t) ~ ^(t), and for 
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Sd galaxies \I/(t) = const, is assumed. Characteristic timescales for SF t* 
(for spirals defined via JjJ* ■ dt = 0.63 • G |t=o) thus range from 1 Gyr for 
ellipticals to 2, 3, 10, and 16 Gyr for Sa, Sb, Sc, and Sd galaxies, respectively. 

The SFHs, together with the IMF, have been chosen as to provide agree- 
ment of our model galaxies after a Hubble time of evolution with integrated 
colours, luminosities, absorption features (E/SOs), emission line strengths (spi- 
rals), and template spectra from the UV through NIR (Kennicutt 1992, Kin- 
ney et al. 1996, cf. Moller et al. 1998), as well as with characteristic HII 
region abundances (:= measured at the effective radius) typical for the re- 
spective galaxy types (Zaritsky et al. 1994, Oey & Kennicutt 1993, van Zee 
et al. 1998, Ferguson et al. 1998). 

2.4 Advantages and Shortcomings 

The combined approach to study the chemical evolution of ISM abundances 
and the spectrophotometric evolution of the stellar population allows - with 
the same number of parameters (IMF and SFH) as any single aspect model 
- for a much larger number of model observables (spectra, luminosities and 
colours from UV through NIR, emission and absorption line strengths, gas 
content and metallicity) to be compared to observations. The combination 
with a cosmological model provides a long redshift or time baseline for com- 
parison with galaxy data, allowing not only to test the models, constrain the 
parameters, but also to make numerous predictions testable by future obser- 
vations (F.-v.A. PhD Thesis 1989). The Gottingen Evolutionary Synthesis 
code has a remarkable analytical potential. It is possible to trace back - at 
any time and in its time or redshift evolution - the luminosity contribution 
to any wavelength band of every single stellar mass, of the various spectral 
types, luminosity classes, and metallicity subpopulations. Ejection rates of 
every individual element, as well, are monitored for every stellar mass, nucle- 
osynthetic origin (PN, SNII, SNIa), and metallicity subpopulation. 

The SFHs we use for different spectral types of galaxies are similar to 
those used by other groups (Bruzual & Chariot 1993, Rocca - Volmerange 
& Guiderdoni 1988) and meant to be global SFRs - averaged over the entire 
galaxy and reasonable intervals of time. SFRs fluctuating around our smooth 
SFHs - either locally or on short timescales - could not be discriminated from 
their smooth idealisations in integrated galaxy properties after sufficiently 
long lookback times. 

Our models are simple 1-zone descriptions without any dynamics or spatial 
resolution, meant to describe global average quantities like integrated spectra, 
luminosities, colours, emission and absorption line strengths, or characteristic 
HII region abundances as measured around the effective radius. While the 
finite lifetimes of individual stars before they restore their partly enriched 
material to the ISM are properly accounted for in the chemical modelling (we 
do NOT use an Instantaneous Recycling Approximation), the mixing of the 
gas is assumed to happen instantaneously. If not indicated otherwise, models 
are closed boxes without any inflow (or outflow). While, clearly, real galax- 
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ies may not be closed boxes over cosmological timescales, this simplification 
reflects our ignorance of the time or redshift evolution of gas infall rates and 
infall abundances. Those cannot yet be taken from hierarchical galaxy forma- 
tion models that are restricted to Dark Matter. The only consistent way to 
account for mass (and energy) exchange beween a (proto-)galaxy and its en- 
vironment seems to be a coupling of our chemo-spectrophotometric evolution 
with models for cosmological structure and galaxy formation (see Contardo, 
F.-v. A. & Steinmetz 1998 for a first attempt). 



2.5 Results for Nearby Galaxies 

Our cc spectrophotometric evolution models were first presented in Moller et 
al. 1997 where we discuss the comparison between the (mass-weighted) ISM 
metallicities for various galaxy types and the luminosity-weighted metallic- 
ities of their stellar populations as seen in different wavelength bands. At 
late stages, stars in models with const. SFR (Sd) show a stellar metallicity 
distribution strongly peaked at |Z at all wavelengths and close to the ISM 
metallicity. Elliptical models, in agreement with observations, show broad 
stellar metallicity distributions extending from Z = 10 -4 to Z = 0.05 in all 
bands (cf. Fig. 1 in F.-v. A. 1999b). 

In collaboration with D. Calzetti (STScI) we are currently working on a 
consistent inclusion of dust, tying the amount of dust to the evolution of both 
the gas content and the metallicity in our models, as well as including different 
spatial distributions of dust and stars in different galaxy types (cf. Moller et 



al. astro-ph / 9906328 for first results). 



3 CC Spectro - Cosmological Evolution: 
Selected Results 

Before the spectrophotometric evolution as a function of redshift can be stud- 
ied for any given set of cosmological parameters, evolutionary corrections e\ 
(since a galaxy at z > is younger) and cosmological corrections - also called 
k-corrections - k\ (since its spectrum is redshifted) are required. 

For redshifts z 3> 1 the attenuation of galaxy light by the cumulative effect 
of intervening hydrogen has to be taken into account. Intergalactic neutral 
hydrogen HI largely comes in the form of Lya clouds causing a forest of narrow 
low column density absorption lines in the featureless continua of background 
QSOs. The cumulative effect of intergalactic HI distributed stochastically 
along the line of sight to galaxies at z > 2.5 has been shown by Madau (1995) 
to significantly attenuate the spectra of distant galaxies at restframe wave- 
lengths below 1216 A. The average attenuation obtained from a large sample 
of sightlines given by Madau et al. 1996 is applied to our redshifted model 
galaxy spectra. It additionally weakens their fluxes below restframe A = 1216 
A and its effect is included in our cosmological corrections. 

The age of a galaxy at redshift z = is given by 
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to := t g al(z — 0) :— tHubblc(z — 0) — tHubble(zf) 



The red colours of present-day elliptical galaxies require galaxy ages in the 
range 12 - 15 Gyr. Reasonable combinations of the cosmological parameters 
are therefore restricted by this minimum age requirement. 

Model galaxy spectra convolved with filter response functions yield abso- 
lute fluxes and magnitudes. Luminosities at z = are normalised to the av- 
erage absolute B - band luminosities of the respective galaxy types observed 
in Virgo (cf. Sandage et al. 1985) before redshifting the synthetic spectra. 

Apparent magnitudes m A of high redshift galaxies in any filter A are ob- 
tained from model absolute magnitudes M A via the bolometric distance mod- 
ulus BDM(H , o , A ), evolutionary and cosmological corrections: 

m A (z) = M A (z = 0, t ) + BDM(z) + e A (z) + k A (z) 

The cosmological or k - correction k A in any wavelength band A accounts for 
the magnitude difference between a galaxy of age t locally and the same 
galaxy spectrum redshifted to z 

k A (z) :=M A (z, t ) -M A (0, t ) 

and can also be calculated for observed galaxy spectra. In this case, the 
maximum redshift to which this is possible depends on how far into the UV 
the observed spectrum extends. Our model galaxy spectra at z = extend 
from 90 A through 160 and, hence, allow for cosmological corrections in 
optical bands up to z > 10. Evolutionary corrections e A account for the age 
difference between a galaxy of today's age t Q and the same galaxy at a younger 
age tg a i (z) corresponding to the redshift z when its light was emitted 

e A (z) := M A (z,t ga i(z)) - M A (z, t c ) 

Evolutionary corrections, of course, cannot be given without an evolutionary 
synthesis model. It is important to stress that both the evolutionary and the 
cosmological corrections do not only depend on the cosmological parameters 
but also on the SFH or spectral type of the galaxy. 

In Moller et al. 2000a (in prep.) the redshift evolution of cc galaxy spectra, 
evolutionary and cosmological corrections, apparent magnitudes, and colours 
will be presented for two different cosmologies. 

3.1 CC Models Compared to Solar Metallicity Models 

A first comparison of the spectrophotometric evolution of model galaxies de- 
scribed in a cc way with those described using only solar metallicity input 
physics is presented in Moller et al. 1997. 

The most important result is that cc spiral models, as compared to Z Q 
models, are brighter by 1 - 2 mag in B and < 1.5 mag in K at redshifts z > 1. 

This has significant implications for the interpretation of high redshift 
galaxy data, e.g. for our understanding of the Faint Blue Galaxy Excess. 
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It means that - contrary to earlier expectations - normal spiral galaxies at 
z ~ 1 . . . > 2 can still contribute to galaxy counts around B ~ 27 mag, and if 
counts are compared to cc models, less of an excess is expected. The impact of 
cc modelling on the Luminosity Functions (LFs) in various wavelength bands 
and on the redshift evolution of LFs of specific galaxy types is currently being 
explored. 

3.2 Comparing with Lyman Break Galaxies 

In any type of galaxy at any time, the flux is close to zero below the Lyman 
break at 912 A. While galaxies with passively evolving stellar populations do 
not have important UV fluxes longward of 912 A either, those with active SF 
do show significant UV fluxes longward, while being self-absorbed shortward 
of 912 A. At z > 2.5 the Lyman break has moved beyond the U - band, 
causing star forming galaxies to drop out of deep U images while visible in 
B, V, etc. images. By z > 3.5 the Lyman break has moved beyond the V - 
band, causing star forming galaxies at z > 3.5 to also drop out of the V - band 
while visible at longer wavelengths. Per deGnitionem only galaxies with active 
SF are detected with this technique. A passive galaxy at redshift z = 3, if it 
existed, would neither be detected in U nor in V or any other optical band. 
Its restframe flux steeply increasing around 4000 A only, it would first appear 
in the NIR H-band. 

Specific colour criteria have been developed to isolate galaxy candidates 
at 2.5 < z < 3.5 and at 3.5 < z < 4.5, called Lyman Break Galaxies (LBGs) 
(see e.g. Steidel et al. 1995, 1999). Follow-up spectroscopy has proven the 
efficiency of this technique and provided fairly large samples of spectroscopi- 
cally confirmed galaxies by today: > 560 galaxies at z ~ 3 and > 46 galaxies 
at z ~ 4 both from the Hubble Deep Field and ground based deep imaging 
surveys. The FORS Deep Field currently observed at the ESO-VLT will soon 
increase the number of very high redshift galaxies by a significant factor. 

The nature of the LBGs is controversial. Appearing quite compact with 
scale radii of 1 - 3 kpc they were suspected to be the progenitors of present- 
day spheroidal galaxies or bulges (Steidel et al. 1996, Giavalisco et al. 1996, 
Friaga & Terlevich 1999) or subgalactic fragments (Somerville et al. 1998). 
Caution seems required, however, comparing rest frame UV scale lengths of 
LBGs with optical scale lengths of local galaxies, since local galaxies on UV 
images look very different in overall morphology - not to mention scale lengths 
- from how they look in the optical (cf. Hibbard & Vacca 1997). The high 
surface density of LBGs, their SFRs - estimated from UV fluxes to be in the 
range 3-60 M© yr _1 -, and luminosities are easier to understand if they are 
the progenitors of spirals. 

Comparing the redshift evolution of luminosities and colours of our cc 
galaxy models with the first set of spectroscopically confirmed HDF dropout 
galaxies (Lowenthal et al. 1997) in Fig. 1, we find that all these LBGs at 
2 z < 3.5 are well compatible with normal spiral galaxy progenitors. 

SFRs derived for LBGs from KECK spectroscopy - corrected for moderate 
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Figure 1: Luminosity and colour evolution of cc spiral models with and with- 
out attenuation compared to spectroscopically confirmed LBGs. 

dust extinction and consistent with their non-detection with SCUBA (Chap- 
man et al. 1999) - also agree well with spiral model SFRs in this redshift 
range. These same results are obtained for two different cosmological mod- 
els (H = 50, q Q = 0.5) and (H Q = 75, q Q = 0.05), with z f = 5. The redshift 
of galaxy formation does not have any visible effect on model luminosities 
or colours in this redshift range unless Zf < 5 for these galaxies. If formed 
around z <~ 5, classical initial collapse elliptical galaxy models are significantly 
brighter and have SFRs > 20 times higher than LBGs. An extensive compar- 
ison including all presently available LBGs is in preparation. 

Clearly, our assumption of one single value for the redshift of galaxy for- 
mation is an oversimplification as cosmological structure formation models 
rather show protracted epochs, extending to very low redshifts for the forma- 
tion of some galaxy types or masses. No generally accepted specific prediction 
for the visible parts of galaxies, however, seems to have emerged so far. 

4 SSPs for Various Metallicities 

The time evolution of SSPs of different metallicities, i.e. of stellar systems 
formed in one short burst of star formation (tsF = 10 5 yr) with one age and 
one metallicity, is useful not only for the interpretation of star clusters, but 
also for combination with various kinds of dynamical galaxy evolution or 
cosmological galaxy formation models, once a SF criterium is specified. Any 
kind of extended SFH in a galaxy or some part of it - as complicated as it 
might be - is readily expanded into a series of single bursts or SSPs. 

In SSPs (as opposed to galaxies with smooth SFHs), the discreteness of 
the stellar mass spectrum with evolutionary tracks available causes strong 
fluctuations in luminosities and colours, that require a posteriori smoothing. 
To avoid this without interpolating tracks or using isochrones we developed a 
Monte Carlo Approach. 
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In principle, this method is simple. The HRD population at any timestep is 
obtained by chosing at random a large number of stellar masses, interpolating 
their lifetimes, and using for each mass not contained in the original track table 
appropriate proportions of the two adjacent tracks with each time interval 
along those two tracks increased or decreased by factors calculated such that 
the lifetime of the "artificial star" equals the interpolated lifetime. 

In Kurth et al. 1999, we present results from our Monte Carlo SSPs. Evo- 
lution of luminosities UBVRIJHK, colours, and absorption indices is shown 
for SSPs of 6 different metallicities —2.3 < [Fe/H] < +0.5 over age ranges 
from 10 7 yr to 16 Gyr. While early in the evolution of a star cluster, changes 
in broad band luminosities and colours are generally rapid, they get weaker 
and weaker with increasing age. SSP models are compared to observations of 
globular clusters in the Milky Way and M32 and theoretical calibrations for 
various indices in ter ms of [Fe/H] are presented in their time evolution. Re- 
sults are available at http:/ / www, uni-sw.gwdg. de/~ ufritze/okurth / ssp.htm, . 

Work on the detailed spectral evolution of SSPs from the UV through the 
NIR is currently in progress. 



5 Chemically Consistent Chemical Evolution 

Timmes et al. 1995 and Portinari et al. 1998 used stellar yields for a range of 
metallicities to model the chemical evolution of the Milky Way and the solar 
neighbourhood. In the following, we will present some of the results from a 
comparison of our cc chemical evolution models with observed abundances in 
the ISM of Damped Lya Absorbers (DLAs). 

5.1 CC Spiral Models compared to DLAs 

DLAs show radiation damped Lya absorption lines in the spectra of back- 
ground QSOs. These damped Lya lines are due to high column density gas 
(log N(HI) [cm -2 ] > 20.3) and usually accompanied by a large number of low 
ionisation lines of Al, Si, S, Cr, Mn, Fe, Ni, Zn, ... with the same absorption 
redshift. Only if the often complex velocity structure in the lines can be fully 
resolved, as e.g. in high-resolution spectroscopy on KECK and WHT with 
AA/A = 60 000, precise element abundances can be derived. These are be- 
coming available for a large number of DLAs over the redshift range ... > 4 
(Boisse et al. 1998, Lu et al. 1993, 1996, Pettini et al. 1994, 1999, Prochaska 

6 Wolfe 1997, ...). 

Based on similarities between their HI column densities and those of local 
spiral disks, between their comoving gas densities at high z and densities of gas 
+ stars in local galaxies, and based on line asymmetries indicative of rotation, 
DLAs are though to be (proto-)galactic disks along the line of sight to distant 
QSOs (e.g. Wolfe 1995, Prochaska & Wolfe 1997, 1998, Wolfe & Prochaska 
1998). Alternatively, Matteucci et al. 1997 propose starbursting dwarf galax- 
ies, Jimenez et al. 1999 LSB galaxies, and Haehnelt et al. 1998 subgalactic 
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Figure 2: Redshift evolution of [Zn/H] (2a) and [Fe/H] (2b) for models Sa 
and Sd compared to observed DLA abundances. The straight line in Fig. 2b 
shows where the model gas content drops below 50%. 

fragments to explain DLA galaxies at low and high redshifts, respectively. 

For any cosmological model the time evolution of metallicity Z, element 
abundances [Xi/H], gas content G, and SFR directly transforms into the cor- 
responding redshift evolution. We stress that for a given SFH and IMF our 
models yield absolute abundances that do not require any scaling or normali- 
sation. After referring all observed DLA abundances to one homogeneous set 
of oscillator strengths and solar reference values, we compare them with our 
cc chemo-cosmological models in Lindner et al. 1999. 

For a standard cosmological model (H Q = 50, o = 1, A Q = 0, Zf = 5) it 
is seen on the example of Zn, which locally is known not to be depleted on 
dust grains, that our Sa and Sd models bracket the redshift evolution of DLA 
abundances from z > 4.4 to z ~ 0.4. The entirety of Zn abundances observed 
in DLAs fall between an upper envelope provided by the rapidly enriching Sa 
model and a lower envelope made up by the slowly enriching Sd model. The 
enrichment evolution of Sb and Sc models run between those for Sa and Sd 
(cf Fig. 2a) . Similar agreement is found for all 8 elements with a reasonable 
number of DLA abundances, i.e. for Zn, Fe, Si, Cr, Ni, S, Al, Mn. Phillipps & 
Edmunds 1996 and Edmunds & Phillipps 1997 argue that the probability for 
an arbitrary QSO sightline to cut through an intervening gas disk and produce 
DLA absorption is highest around the effective radius where our models were 
shown to match average HII abundances in the respective spiral types by 
z = 0. Hence, our models bridge the gap from high-z DLAs to nearby spirals. 
We conclude that from the point of view of abundance evolution, 
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DLA galaxies may well be the progenitors of normal spirals Sa — Sd, 

although we cannot exclude that some starbursting dwarf or LSB galaxies may 
also be among the DLA galaxy sample. If, at the highest redshifts, damped 
Lya absorption were caused by subgalactic fragments still bound to merge, our 
simplified 1-zone models are meant to describe the total SFR of all fragments 
ending up in one galaxy by low redshift. The weak redshift evolution of DLA 
abundances is a natural result of the long SF timescales for disk galaxies. 
The range of SF timescales t* for spirals from Sa through Sd fully explains 
the abundance scatter observed among DLAs at any redshift. 

Note that our comparison of spiral galaxy models with DLA abundance 
data extends to redshifts z > 4.4, i.e. over lookback times of > 90% of the 
age of the Universe. 

In comparison with solar metallicity models the influence of the metallic- 
ity dependent yields is seen to vary from element to element. Whenever a 
significant difference is seen, the cc models give better agreement with the 
data than Z Q models. Changing the SN explosion energy (yields for model C 
from Woosley & Weaver) has a minor effects (see also Fig. 7a-h in Lindner et 
al. 1999). 

Somewhat surprisingly, elements which locally deplete strongly onto dust 
grains (like Fe or Cr) are as well described by our models as are non-refractory 
elements like Zn (cf. Fig. 2b). We hesitate, however, to draw conclusions 
about the importance of dust in DLAs in view of the uncertainties in the 
stellar yields and the simplicity of our closed-box models. 

We also tried models with constant infall rates and primordial infall abun- 
dances. In this case, SFHs have to be adjusted as to still give agreement at 
z = with average type-specific HII region abundances, colours, etc. This is 
the reason for moderate constant infall rates to only marginally change the 
redshift evolution of model abundances and to not affect any of our conclu- 
sions. A redshift dependent infall rate, however, with or without evolving 
infall abundance, might affect our results. The difficulty is to constrain the 
additional parameter (s) without embedding the galaxy into a cosmological 
context. Mass estimates from column densities, linear dimensions (DLAs 
seen in close QSO pairs), and rotation velocities (Prochaska & Wolfe 1997) 
indicate that - at least some of - these systems at z ~ 2 — 3 already have the 
masses of local spirals and, hence, may not require important infall. 

5.2 A Change with Redshift in the DLA Population ? 

While at high redshift all spiral types Sa through Sd seem to give rise to DLA 
absorption, no more data points at z < 1.5 reach close to our early type spiral 
models. The boundary to the observed abundances coincides with the 50 % 
gas-to-total-mass ratio in our models. I.e., by low redshift, the gas poor early 
type spirals seem to drop out of DLA samples. A deficiency of high N(HI) 
DLA systems at low z had been noted before (Lanzetta et al. 1997) and 
attributed to their high metallicity and dust content (Steidel et al. 1997). We 
give an additional reason: as the global gas content drops, the probability for 
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an arbitrary line of sight to a QSO to cut through a high column density part 
of the galaxy decreases, i.e. the cross section for damped Lya absorption 
gets reduced. 

5.3 Spectroscopic Predictions for DLA Galaxies 
and Future Prospects 

Despite considerable efforts and large amounts of telescope time devoted to 
the optical identification of galaxies responsible for damped Lya absorption, 
only less than a handful have been found by today. Contrary to earlier ex- 
pectations, the success rate has not increased when low redshift DLAs were 
found in UV spectra of QSOs. If, indeed, there were a shift in the DLA galaxy 
population towards later spiral types at lower redshift, then this is what our 
models predict since, locally on average, Sd galaxies are fainter by ~ 2 mag in 
B than Sa's. Hence, the low-z late type DLA galaxies should be about as faint 
in B, 1Z, and K as the brightest members (i.e. the early spiral types) of the 
high-z population. Both, an average Sd model at z ~ 0.5 and an average Sa 
model at z ~ 2 — 3 have B ~ 25, 7\L~ 24.5, K <~ 22 mag. Luminosities of the 
few optically identified DLA galaxies (and candidates) are in good agreement 
with our predictions. In particular, if there were early type spirals among 
the low redshift DLA galaxy population they would have B ~ 22.5, 21, 
K ~ 18.5 mag at z~ 0.5 and they would have been detected. Their non- 
detection is consistent with our finding of a change with redshift in the DLA 
galaxy population (see F.-v.A. et al. 1999a, b, c for details). 

DLA galaxies are within the reach of lOm-class telescopes up to redshifts 
z > 3 and trace the normal galaxy population to these high redshifts with- 
out any bias as to high luminosity, radio power, or the like. With information 
about ISM abundances from the metal absorption lines and spectrophotomet- 
ry properties of the stellar population, they will allow for powerful constraints 
on model parameters, ages, and cosmological parameters (cf. Lindner et al. 
1996). Accurate abundance data in very low metallicity DLAs may provide 
clues for the nucleosynthesis at low metallicity. 

The knowledge about the abundance evolution of gaseous spiral disks as 
a function of redshift will allow to predict abundances of stars, star clusters, 
and Tidal Dwarf Galaxies that form in the powerful starbursts accompanying 
spiral galaxy mergers in the local universe and in the past (cf. Sect. 6.3 and 
6.4). As will be shown, knowing these abundances is important for properly 
interpreting observed colours and luminosities of star clusters and Tidal Dwarf 
Galaxies in terms of ages, masses, etc. Spectroscopy of those objects, in turn, 
provides an independent cross check of our abundance predictions. 



6 Interactions and Starbursts 
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6.1 Importance of Galaxy Interactions and Starbursts 

Evidence for the importance of galaxy interactions and merging is coming 
both from theoretical galaxy formation modelling and observations. 

In a variety of cosmological contexts (CDM with and without A, mixed 
DM), semianalytic as well as numerical structure formation models predict 
galaxies to build up hierarchically from a series of mergers of successively 
larger masses (see e.g. Carlberg 1990, Lacey & Cole 1993, Kauffmann & 
White 1993). With increasing redshift, observations indicate increasing num- 
bers of physically close galaxy pairs (e.g. Zepf & Koo 1989), galaxies seen in 
interaction, and disturbed galaxy morphologies indicative of a recent inter- 
action (e.g. Conselice & Bershady 1999). All strong galaxy interactions or 
mergers are observed to be accompanied by bursts of star formation if one or 
both of the galaxies are gas rich. When falling into a cluster, the spiral rich 
field galaxy population must somehow be transformed into the E/SO/dE rich 
cluster population. These transformation processes also involve interactions 
- among individual galaxies or of an infalling galaxy with the cluster envi- 
ronment, potential, or its hot intracluster medium. In gas rich galaxies some 
of these processes, as well, may trigger starbursts - sometimes powerful and 
global, as reported for E+A galaxies by Poggianti & Barbaro 1996. 

In dissipationless stellardynamical mergers (of disks or spheroidal sys- 
tems), the remnant after partially complete violent relaxation usually is a 
spherical system. Its central phase space density cannot significantly increase 
beyond that of the progenitors. Therefore, it does not seem possible to pro- 
duce the high central density cores of massive ellipticals by dissipationless 
mergers of less massive subcomponents or dwarf galaxies which intrinsically 
have lower central densities. Gaseous mergers, however, are highly dissipative, 
the gas is driven to the center very efficiently, and gas densities comparable to 
the central stellar densities in massive ellipticals are reached in high-resolution 
simulations. In Ultraluminous Infrared Galaxies (ULIRGs), all of which have 
been shown to be mergers of gas rich spirals, central molecular gas densities 
of the order of 1000 M Q pc -3 are indeed observed. 

It has not been possible yet to treat the full dynamics of mergers involv- 
ing stars, DM, and gas in at least 3 phases (Hot X gas, HI, cold molecular 
gas), and to consistently include SF and feedback (under extreme conditions) 
into high-resolution stellar + gaseous dynamical models of interacting gas rich 
galaxies. It is clear, however, that strong bursts of SF or/and efficient fuelling 
of a central AGN may occur. Which alternative is chosen or dominates may 
depend on the properties of the galaxies involved, the geometry of the en- 
counter, and possibly on the stage of the interaction (see e.g. F.-v.A. 1996 for 
a review in Galaxy Interactions). 

By studying interaction-triggered starbursts with spectrophotometric and 
chemical evolutionary synthesis models we hope to learn about the SF process, 
SF efficiencies, etc. under the violent conditions in mergers/strong interac- 
tions that seem to significantly differ from the comparatively "peaceful" SF 
situation in our Milky Way. 
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6.2 Selected Results 



We calculated an extensive grid of starburst models. In the idealisation of 
equal-type mergers Sa-Sa, Sb-Sb, Sc-Sc, Sd-Sd, occuring at a variety of evolu- 
tionary stages of the parent galaxies, we model starbursts of various strengths 
and durations on top of the pre-existing stellar population. The age and type 
of the progenitor galaxies determine the luminosities and colours of the stellar 
population and the ISM metallicity at the onset of the burst. The gas content 
of the spirals at merging sets an upper limit to the gas reservoir available for 
SF in the burst. Burst strengths being defined as the increase of the stellar 
mass S during the burst, b := o ASbur3t , this means that e.g. Sa-Sa mergers 

•^pre — burst 

in the local universe cannot have bursts as strong as Sd-Sd mergers today or 
Sa-Sa mergers earlier in their evolution when the galaxies contained more gas. 

The evolution of galaxies before, during, and after the burst was studied 
in terms of two-colour diagrams, as well as in their luminosity and metallicity 
evolution for solar metallicity models in F.-v.A. & Gerhard 1994a. 

The models show that after a burst spiral-spiral mergers of galaxy types Sa, 
Sb, Sc occuring 3-4 Gyr ago readily develop typical elliptical galaxy colours 
by today - provided SF ceases completely. Sd mergers remain too blue for 
> 4 Gyr. In detail, the timescale until a merger remnant with its postburst 
reaches the observed colour range of E/SOs depends on the progenitor types, 
the burst strength, and on the wavelength range considered. 

Using our grid of models we analysed the starburst in the prototypical 
merger remnent NGC 7252, for which a wealth of additional information is 
available (e.g. spectra from UV - optical, HI-, HII-, CO-maps, dynamical 
modelling). Despite its enormous tidal tails it already features an azimutally 
averaged r 1 / 4 - profile in its inner part, and it is the oldest from Toomre's 
(1977) dynamical age sequence of interacting galaxies. Comparing broad band 
colours UBVR with our model grid, we were left with some range of possible 
combinations of burst strengths and ages. Within this cell of parameter space, 
we in detail compared our model spectra with the spectrum of NGC 7252 
showing deep Balmer absorption lines and a small emission component. 
This allowed us to identify a very strong starburst that increased the stellar 
population by < 50 % and started ~ 1 Gyr ago, i.e. around the time of 
1 st pericenter as obtained from dynamical modelling (cf. Hibbard & Mihos 
1995). Both in terms of morphology and colours, NGC 7252 could evolve into 
an elliptical galaxy within the next 1-3 Gyr, provided its SFR — > 0. Its 
present centrally concentrated SFR as obtained from an IUE spectrum (F.- 
v.A. et al. in prep.) is ~ 3 M Q yr~ 1 , consistent with the weak emission. 
HI falling back from the tidal tails (Hibbard & Mihos 1995) together with gas 
restored by dying burst stars supports the present SFR and may do so for 
another few Gyrs, allowing NGC 7252 to rebuild itself a small disk of stars 
and gas. In this case, both in terms of morphology and colours, it might 
rather come to resemble an SO or even an Sa galaxy. 

In an attempt to identify typical values and ranges for burst parameters 
and study their possible dependence on galaxy and interaction properties, 
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Liu & Kennicutt's (1995) sample of interacting galaxies was analysed by 0. 
Kurth in his diploma thesis (Gottingen 1996). Only on the basis of those 
typical values and relations can interaction induced starbursts consistently be 
included into the statistical evolution of high redshift galaxies, as e.g. into 
models for the redshift evolution of the luminosity function or into cosmolog- 
ical galaxy formation scenarios. In general, burst durations can only weakly 
be constrained to be of the order of a few 10 8 yr for spiral - spiral mergers, i.e. 
definitely longer than for Blue Compact Dwarf Galaxies (BCDs), where they 
are of the order of 10 6 yr. In both cases, burst durations reflect the dynamical 
timescales of the systems involved. One of the most important results is that 
starbursts in galaxy mergers can be very strong, increasing the stellar 
mass by as much as 30 % or even more. In all massive interacting systems, 
burst strengths > 10% were found, larger by a factor 10 — 100 than the burst 
strengths obtained for BCDs (Kriiger et al. 1995). 

For NGC 7252, the peak SFR during the burst must have been several 
hundred M Q yr _1 , comparable to those of ULIRGs. Assuming a gas content 
from the high end of what is typical for Sc galaxies, the SF efficiency (SFE) 
- defined as the mass of stars formed in the burst relative to the mass of gas 
available r\ := G AS|J ° rat - - must have been very high, of order 30 - 45 %, a 
factor > 10 larger than in the Milky Way (F.-v.A. & Gerhard 1994b). 

CO, HCN, and CS observations of ULIRGs tracing molecular gas at den- 
sities n ~ 500, 10 , 10 5 cm~ 3 , respectively, explain their high SFEs as a con- 
sequence of their high fractions of molecular gas at the highest densities 
(Solomon et al. 1992). While for BCDs a trend for burst strengths to decrease 
with increasing galaxy mass was found (Kriiger et al. 1995), burst strengths 
and SFEs in interacting galaxies with masses 1-2 orders of magnitude higher 
are generally higher by the same 1-2 orders of magnitude. Does this mean, 
that the SF process is intrinsically different in different environments ? (See 
Sect. 5 in F.-v.A. 1996 for a discussion of "violent vs peaceful" SF and the 
possible relation with differences in the molecular cloud structure). 

SFEs that high are usually assumed for the collapse of protogalactic clouds 
in the early Universe that gives rise to globular cluster formation. We specu- 
lated that globular cluster formation might have been possible in the starburst 
of NGC 7252 and predicted the metallicity of both stars and clusters forming 
in the burst to be Z* > Z^J 1 Gyr ^ > | • Z Q . 

6.3 Young Star Clusters in Mergers 

Within a few years, bright Young Star Cluster (YSC) populations were de- 
tected with HST in NGC 7252 and a number of other interacting galaxies and 
merger remnants. They are interesting for many reasons. E.g., star clusters 
are better suited than the integrated light to study the age and duration of 
the starburst. Being most probably observed around the same galactocentric 
radii at which they formed, they thus allow to study the spatial extent of the 
starburst and - for sufficiently large cluster populations - its time evolution. 
YSCs offer the possibility to study star and cluster formation processes in 
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a violent environment, in particular in comparison with observations of the 
molecular gas content or cloud structure. Using SSP models for the inter- 
pretation of YSC systems in the old merger remnant NGC 7252 (F.-v.A. & 
Burkert 1995) and in the young interacting Antennae galaxies NGC 4038/39 
(F.-v.A. 1998a), we show that metallicities of YSCs are crucial to age-date 
them on the basis of observed broad band colours and to predict their fu- 
ture luminosity and colour evolution. For a few of the brightest YSCs in 
NGC 7252, our metallicity prediction from spiral models was confirmed by 
spectroscopy (Schweizer & Seitzer 1993, 1998). 

The question if some, many, or most of these YSCs could be young Globular 
Clusters (GCs) is of eminent importance. If a secondary GC population, 
comparable in number to the original one, can be formed in a merger induced 
starburst, this would invalidate the last surviving argument against the for- 
mation of - at least some - elliptical galaxies from spiral-spiral mergers based 
on the difference in specific GC frequencies between elliptical galaxies and 
spirals. Analysing WFPC1 data for some 40 YSCs in NGC 7252 and > 500 in 
NGC 4038/39 we tentatively conclude that the bulk of the YSC populations 
in both systems may well evolve into "normal" GC populations in terms of 
colour distributions, luminosity and mass functions (F.-v.A. 1999a). Before 
definite conclusions can be reached, we are currently reanalysing WFPC2 data 
and studying YSCs in an age sequence of interacting, merging and merged 
galaxies to assess the impact of dynamical cluster destruction processes. 

Because of its enhanced metallicity, a secondary GC population might 
eternally testify to a merger origin, still at times when tidal tails, kinematic 
peculiarities, and fine structure will long have disappeared, and colours will 
no longer reveal a past starburst. Analyses of GC systems in a number of 
elliptical galaxies - from low-luminosity ellipticals all through cD galaxies - 
have revealed bimodal colour distributions in many bright ellipticals, including 
2 SOs, broad or multi-peak colour distributions in all cD galaxies investigated, 
and single-peak distributions in a few low-luminosity ellipticals (Kissler - 
Patig this volume, Kissler - Patig et al. 1998, Gebhardt & Kissler - Patig 
1999). With 10 m telescopes in combination with HST imaging, MOS becomes 
feasible for YSCs and even old GCs out to Virgo distances. In comparison with 
SSP models, it will allow to decompose the colour distribution of star cluster 
systems into metallicity and age distributions and, thereby, give information 
about the formation of their parent galaxies. 

6.4 Tidal Dwarf Galaxies 

Not only secondary populations of YSCs are formed during mergers of gas- 
rich galaxies, but also a new galaxy formation mechanism has been detected 
in these systems a few years ago. In the extended tidal tails of several in- 
teracting systems bright blue star forming knots, often associated with large 
HI concentrations, are observed with masses and luminosities typical of dwarf 
galaxies. Two of these Tidal Dwarf Galaxies (TDGs) are detected in NGC 
7252, one in NGC 4038/39, two in Arp 105, several in the "Superantennae", 
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and other systems (see Due et al. 1998 for a recent review). 

TDGs form from "recycled" material torn out from a spiral disk. They de- 
viate from the dwarf galaxy luminosity - metallicity relation in the sense that 
they have enhanced metallicities for their luminosities as compared to "non- 
recycled" dwarf galaxies. Abundance determinations from their HII region- like 
spectra agree well with predictions from our spiral models. 

While TDGs forming in the lower parts of tidal features are expected to fall 
back into the merger remnant on timescales of 10 8 to few 10 9 yr, those at the 
tips of the tails will probably escape and might survive as independent entities. 
Kinematic independence from the parent galaxy has been demonstrated in a 
few cases (Due et al. 1997, Due & Mirabel 1998). 

Independently, both N-body simulations for the stars and hydrodynami- 
cal models for the gas in merging galaxies feature condensations along their 
respective tidal tails, leading to two competing scenarios for the formation of 
TDGs (cf. Barnes & Hernquist 1992, Elmegreen et al. 1993). 

Evolutionary synthesis models involving fractions of the stellar popula- 
tions of spiral progenitors plus starbursts of various strengths are used to 
analyse the stellar populations and the evolutionary states of TDGs. Gaseous 
continuum and line emission are included. Analysing a small sample of TDGs 
observed by P.-A. Due with a coarse grid of models (F.-v.A. et al. 1998) we 

□young 

find evidence for strong bursts, b :— gold = 0.1 — 0.4, on top of 'old' stellar 
populations, i.e. populations with the age distribution of stars in their parent 
spiral (S young and S old being the masses of young and 'old' stars, respectively). 
Since these strong bursts completely dominate the light in the optical, NIR 
data are required to constrain the mass contribution of the 'old' population. 
The latter is important for the fading that TDGs will experience in the future, 
as well as for their dynamical evolution. Without a significant old population 
the strongest bursts might disrupt the TDG. Both, the fading and future dy- 
namical evolution, in turn, are relevant for a possible cosmological significance 
of TDGs. In the past, the merger rate was higher, galaxies were more gas-rich 
and probably less stable so that the production of a significant population of 
TDGs might be expected. 

To explain that part of the faint blue galaxy excess that is due to dwarf 
galaxies, Babul & Ferguson 1996 invoke a hypothetic population of dwarf 
galaxies, the formation of which is delayed until z ~ 1 by the intergalactic 
radiation field. Properties they require for their dwarf galaxies to explain the 
faint blue galaxy excess are very similar to those of our TDGs. Depending 
on their fading, TDGs might well explain part of the faint blue galaxy excess. 
Moreover, the remnant problem faced by Babul & Ferguson and Ferguson & 
Babul 1998 would be alleviated if part of the TDGs would spiral back into 
the merger remnant on timescales of a 10 s — 10 9 yr. 

In Wcilbacher et al. 2000 (submitted) we analyse a sample of 10 inter- 
acting galaxies from the AM catalogue with imaging in B, V, R. Comparing 
colours of knots along tidal structures with an extensive grid of TDG model 
calculations, a number of promising TDG candidates are identified. Knots 
with colours not explained by models most probably are background galaxies. 
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TDG candidates in this first sample have young burst ages of ~ 7 Myr, burst 
strengths b < 0.2, and will fade by up to 2.5 mag in B, on average, within 200 
Myr after the burst. Follow-up spectroscopy with VLT and HET is underway, 
models are currently being calculated for their detailed spectral evolution. In 
any case, this recently discovered mode of ongoing galaxy formation in the 
local Universe from recycled gas and stars is an interesting field of study just 
at the beginning of being explored. 

7 Conclusion and Outlook 

I presented a versatile model for the evolution of stellar populations and gas 
that offers a variety of applications from star clusters to nearby and distant 
galaxies. Only a few of them have been presented here. With a minimum 
number of parameters our combined chemical, spectrophotometric and cos- 
mological evolution models describe a large number of observables and provide 
a long evolutionary baseline to compare with high-redshift galaxy abundances 
and spectra and understand the evolution of various galaxy types from the 
earliest phases to the present. The chemically consistent treatment is a first 
attempt to consistently combine 2 out of 3 aspects of galaxy evolution that 
nature, too, has coupled. 

Interactions play an important role for the evolution of galaxies over cos- 
mological timescales. While in its present state, the model does not include 
any dynamical aspect nor spatial resolution, we tried and studied the effects 
of starbursts accompanying galaxy interactions if gas is involved. Several in- 
teresting phenomena were observed in this context, as e.g. the formation of 
large populations of bright star clusters and of "recycled" Tidal Dwarf Galax- 
ies. Application of our models provided a first step to understand the nature 
and properties of these systems as well as their possible future evolution. 

The model has allowed for a series of precise observational predictions, 
part of which became verified already while others keep standing for a test. 

Our first attempt to also include the 3 rd aspect of galaxy evolution, the 
formation and dynamical evolution of galaxies in their cosmological environ- 
ments, could not be discussed here (cf. Contardo et al. 1998). 

Over the years this model was developed, extended and refined, with its 
input physics continuously updated, observational extragalactic research has 
seen a tremendous progress. The amount of information from HST and large 
ground-based surveys is enormous and several quantum steps have been per- 
formed, e.g. concerning image resolution with HST, spectral resolution with 
KECK and WHT, and the number of high-redshift galaxies by the Lyman 
break technique. With 10 m telescopes, like VLT and HET, observational 
progress is challenging theory to keep path. The particularly close interplay 
between observations and the conceptually simple galaxy evolution modelling 
presented here, has proven very fruitful and stimulating for both sides. 
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